Bragg spectroscopy of trapped one dimensional strongly interacting bosons in optical lattices: 

Probing the cake-structure 



Guido Pupillo^, Ana Maria Rey^ and Ghassan George Batrouni'^ 
^ Institute for Quantum Optics and Quantum Information of the Austrian Academy of Sciences, 6020 Innsbruck, Austria 
^Institute for Theoretical Atomic, Molecular and Optical Physics, 
Han'ard-Smithsonian Center of Astrophysics, Cambridge, MA, 02138, USA and 
^Institut Non-Lineaire de Nice, Universite de Nice-Sophia Antipolis, 1361 Route des Lucioles, 06560 Valbonne, France 

(Dated: February 6, 2008) 

We study Bragg spectroscopy of strongly interacting one dimensional bosons loaded in an optical lattice plus 
an additional parabolic potential. We calculate the dynamic structure factor by using Monte Carlo simulations 
for the Bose-Hubbard Hamiltonian, exact diagonalizations and the results of a recently introduced effective 
fermionization (EF) model. We find that, due to the system's inhomogeneity, the excitation spectrum exhibits 
a multi-branched structure, whose origin is related to the presence of superfluid regions with different densities 
in the atomic distribution. We thus suggest that Bragg spectroscopy in the linear regime can be used as an 
experimental tool to unveil the shell structure of alternating Mott insulator and superfluid phases characteristic 
of trapped bosons. 
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I. INTRODUCTION 



Cold atoms in optical lattices provide a way for realizing 
interacting many body systems in essentially defect free lat- 
tices Paradigms of strongly correlated phenomena such as 
the superfluid (SF) to Mott insulator (MI) quantum phase tran- 
sition have been realized in a three-dimensional lattice thanks 
to the successful application of atom optics techniques to tra- 
ditional condensed matter systems |2]|. Recently, much inter- 
est has been generated by experiments in reduced dimension- 
ahty iHBiSf^i,^. where the role of interactions and quan- 
tum fluctuations is enhanced. Of particular relevance have 
been the realization of a SF/MI quantum phase transition in 
one-dimension (lD)JHJp, and the observation of a gas of 
hard core bosons yl or Tonks-Girardeau (TG) gas 0] ■ 

The latter is a characteristic of strongly interacting one dimen- 
sional systems, where the large repulsion between atoms mim- 
ics the Pauli exclusion principle. As a consequence, there is 
a one to one correspondence of the eigenenergies and eigen- 
functions of TG bosons and non-interacting fermions. This 
coiTespondence holds for all local observables, such as the 
atomic density and fluctuations |9|] . 

The presence in experiments of a confining parabolic po- 
tential superimposed on the lattice provides the possibility of 
studying strongly correlated systems at arbitrary densities. In 
fact, the presence of the quadratic potential typically induces 
the coexistence of alternating superfluid and insulating phases 
with on-site lattice densities which can be larger than one 
This creates a shell structure, which is rem- 
iniscent of the MI lobes of the homogeneous phase diagram 
Jisll . In this case, the SF/MI transition is better understood 
as a crossover than as a phase transition flTll . Recently, it 
was shown that a system of strongly interacting ID bosons at 
high densities has the structure of an array of stacked disjoint 
TG gases fl4ll . This decomposition of the system in indepen- 
dent TG gases, or effective fermionization (EF), allows for the 
computation of static properties such as the density profile and 
fluctuations, as well as dynamical properties. In particular, EF 
has been used to explain the microscopic mechanisms respon- 



sible for the decay of the superfluid current of an interacting 
bosonic gas in a periodic potential 0]. 

It remains a challenge to engineer experimental probes 
for atoms in the strongly correlated regime. Information on 
the excitation spectrum has been obtained by using Bragg 
spectroscopy |8], sparking considerable theoretical activity 
Jisl Ha .17. .18J . Unfortunately, the experiment of Ref. 01 
was conducted in a regime far from linear, which prevented 
the direct comparison of the theoretical results with the exper- 
imental data. In particular, in the case of the inhomogeneous 
system it has not been possible to verify the nature of the dou- 
ble peak structure found in Ref. in the low frequency re- 
sponse, by direct comparison to the existing data of Ref. 

In this paper we study Bragg spectroscopy of bosons in 
the periodic plus quadratic potentials in the linear response 
regime. We perform quantum Monte Carlo simulations and 
exact diagonalizations of the Bose-Hubbard Hamiltonian in 
the strongly interacting regime to compute the dynamic struc- 
ture factor. The latter is found to exhibit a multi-branched 
structure at small excitation frequencies [17]. We compare 
these numerical results to the predictions of the EF model, 
finding good agreement. The central result of this paper is 
that this agreement strongly indicates that this multi-branched 
structure at low-frequency is dominated by the excitations of 
the system's superfluid components. That is, the observed 
different excitation branches can be attributed primarily to 
the response of the superfluid components of the various 
layers of the EF model, whose presence is directly linked to 
the existence of a shell structure in the many body density 
profile. When there are at least two superfluid regions with 
different mean particle numbers, as is the case of Ref. (H, this 
result differs substantially from previous results, where the 
presence of the double-peak structure in the response to the 
Bragg perturbation was linked to the presence of particle-hole 
excitations in the Mott phase il7il . Therefore, because of 
the sensitivity to the presence of superfluids with different 
particle numbers, we suggest that Bragg spectroscopy in the 
linear regime is an ideal tool to characterize the shell structure 
of atoms in the quadratic potential (see also lll9ll20ll2lll '). 
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The presentation of the results is organized as follows. In 
Sec. lUwe introduce the Bose-Hubbard Hamiltonian, shortly 
review the physics of strongly interacting bosons in a lattice 
and explain the basic ideas of the effective fermionization 
model. In Sec. |lll]we introduce the dynamic structure fac- 
tor In Sec. llVl we compare the results of the EF model and 
of exact diagonalizations for a small number of bosons in a 
homogenous lattice. The various branches of excitation are 
explained in terms of (low-frequency) excitations of the su- 
perfluid and (higher-frequency) particle-hole excitations. The 
EF model is shown to well reproduce the excitations of the su- 
perfluid. In Sec.|V]we compare the results of quantum Monte- 
Carlo simulations for experimentally realistic systems to the 
results of the EF model. When there are more than one super- 
fluid regions with different particle densities, we show strong 
evidence that the low-frequency excitations are dominated by 
the excitations of the various superfluids. We then propose an 
experiment in the linear response regime which would detect 
the associated shell structure of the many-body density pro- 
file. Finally, the conclusions are presented in Sec. lVIl 



II. BOSE-HUBBARD HAMILTONIAN 

The Bose-Hubbard {BH) Hamiltonian describes N interact- 
ing bosons in a lattice potential of M sites Ql 
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Here aj is the bosonic annihilation operator at site j, and fij = 
OjUj. is the curvature of the confining potential and U and 
J are the on-site interaction and hopping energies. 

In the homogeneous system, = 0, the effective strength 
of the interactions depends on the atomic density N/M, and 
the ratio j ~ U/J. In particular, for N < M the system is 
effectively a strongly interacting TG gas if U is larger than 
the lattice band width 4 J 0]. Thus, 7 = 7^ ~ 4 deter- 
mines the critical interaction strength in low density systems. 
In particular, the SF/MI transition occurs in a homogeneous 
unit filled lattice at 7 ^ 7c. For arbitrarily large densities, 
n — 1 < N/M < n, with n an integer larger than one, the 
system enters the strongly correlated regime when 7 ^ 7cn. 
In these cases standard fermionization techniques are invalid, 
however the low energy physics can still be well reproduced 
by considering the excess M[N — {n — 1)] bosons as non- 
interacting fermions with an effective kinetic energy nJ sit- 
ting on a plateau formed by M (n — 1) atoms frozen in a Mott 
state with exactly ?i — 1 atoms per site. The validity of this 
approximation relies on the fact that states with more than n 
atoms per site are suppressed by a factor on the order of 1 /j. 
In other words, when 7 ^ jcn, the whole many-body system 
can be conveniently visualized as composed of two indepen- 
dent subsystems; the atoms frozen in the Mott state and the 
extra TG bosons 11 411 . 

For most experiments, a parabolic magnetic confining po- 
tential of frequency proportional to fl is also present, so the 
density profile varies across the lattice. In the trapped case the 



conditions for the formation of a Mott state change dramati- 
cally. For example, it is always possible to create a unit-filled 
MI if iV < M, by varying the depth of the lattice or the mag- 
netic trap frequency. In the trivial J = limit, the density 
profile becomes a "cake" structure with maximal occupation 
n* at the trap center. In our model, we view the density as a 
"layer cake" of n stacked horizontal layers. When J = 0, the 
atoms are completely frozen, and each layer of the cake may 
be viewed as an independent Mott state with Nn unit-filled 
sites. 

For moderate values of J > 0, it has been shown theoreti- 
cally that the density profile still has a cake structure of coex- 
isting superfluid and MI phases. In this case, the density pro- 
file can be visualized as being composed of stacked horizontal 
layers, but because atoms are no longer frozen, in general the 
layers are not independent. However, if number fluctuations 
in adjacent horizontal layers do not overlap in space, all lay- 
ers can be treated independently and, as in the homogeneous 
system, standard fermionization techniques can be applied to 
each layer separately. In this situation single-particle solu- 
tions provide expressions for all many-body observables. We 
call this generalization of the Bose-Fermi mapping extended 
fermionization (EF). 

We want to point out that at variance with the original def- 
inition jT4ll . and for the sake of simplicity, we here refer to 
all kinds of strongly interacting gases as EF, independently of 
the density and the presence or absence of an external poten- 
tial. In the low density limit n = 1, EF reduces to standard 
fermionization. 

In the remainder of this paper we focus on this strongly 
correlated regime, where 7 ^ 7cn. Because J decreases ex- 
ponentially with increasing lattice depth, this regime is easily 
attained experimentally. 



in. THE DYNAMIC STRUCTURE FACTOR 

The dynamic structure factor at temperature T is given by: 

5(g,c.) = l^e-^'^'|(z|p,|j)pJ(fi^-i?,+i?,), (2) 

where, Z is the canonical partition function, i3~^ = Tks, 
with ks the Boltzmann constant and T the temperature. Here, 
Ej are the eigenenergies of the many body eigenstates \i), and 

Pq ^ jj J2j e'^'^-' ajaj is the density fluctuation operator with 
hq the quasi-momentum and d the lattice constant. For the 
quantum Monte Carlo simulations we use a World-Une algo- 
rithm at small finite temperature Tks = O.IJ fl7ll23l . while 
for the exact diagonalizations we use standard linear algebra 
techniques, and T can be taken to be zero. Equation Q de- 
scribes the response of the system to Bragg spectroscopy in 
the linear response regime. It implies that the system responds 
whenever the frequency lu of the Bragg perturbation matches 
the energy difference between two eigenstates. 

For 7 ^ 7cn and at low enough energy, eigenstates are 
densely grouped in energy ranges of the order of a few J l24ll . 
separated by an energy of the order of U. Eigenstates in each 
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energy range are linear combinations of Fock states with the 
same number of empty sites, of singly occupied sites, doubly 
occupied sites, and so on. We thus expect that in the strongly 
correlated regime the system's response to the Bragg perturba- 
tion mirrors the "grouped" structure of the many body energy 
spectrum. 

In standard fermionization, S{q, ui) is obtained by comput- 
ing Eq. (|3 for a system of non-interacting fermions i23ll . 



(3) 



Here and are the n*'' single-particle eigenmodes 
and eigenenergies of Eq.Q with hopping energies J respec- 
tively and j is the lattice site index. denotes the 
Fermi-Dirac distribution. 

For the cases when there are more than one atom per site, 
according to the EF model the density profile can be visual- 
ized as being composed of stacked horizontal layers. In the 
parameter regime where number fluctuations in adjacent hori- 
zontal layers do not overlap in space, all layers can be treated 
independently and the dynamical structure factor of the over- 
all system can be approximated by adding the structure fac- 
tors of the n different layers each one with Nn atoms. The 
expression for S{q,uj) in a given n layer is exactly the same 
than Eq.© but replacing the eigenmodes and eigenenergies 
by the single-particle solutions of Eq.Q with hopping ener- 
gies nJ . Also the chemical potential in the Fermi-Dirac dis- 
tributions must be calculated to consistently have an average 
of Nn atoms in each layer. 

The independent addition of the structure factor of each 
layer is justified in the model, because EF takes into account 
only the responses of the superfluids to the Bragg perturba- 
tion. Thus, when the two superfluid regions are spatially well 
separated, matrix elements coupling low-lying excitations of 
the two superfluids in Eq. (|2ji become zero, because of the 
vanishing overlap of the wave functions. 



IV. HOMOGENEOUS LATTICE (!^ = 0) 

In Fig. n we show S{q,uj) as a function of the Bragg fre- 
quency oj, for a homogeneous lattice with M = 6 and periodic 
boundary conditions. The number of atoms is = 5 (panels 
(a) and (c)) and N = 7 (panels (b) and (d)). The temperature 
T is equal to zero, in all the plots. The continuous, dashed 
and dotted lines correspond to different quasi-momenta, with 

= 7r/3, 27r/3 and tt, respectively. Panels (a) and (b) are the 
results of exact diagonalizations with 7 = 10, while panels 
(c) and (d) are the predictions of the EF model for atoms with 
kinetic energy J and 2 J, respectively. In Fig. 1(a) the exact re- 
sults show large peaks for huj/J<4, and smaller excitations 
for 8 < hoj / J < 14. The value 4 J corresponds to the single 
particle band width of a hole, that is of an empty site tunnel- 
ing in the lattice. These low frequency excitations are thus due 
to the coupling between the ground state and eigenstates with 



at most one atom per site. The smaller peaks in the interval 
8 < ftw/ J < 14 are instead due to the coupling to particle- 
hole excitations (1-ph), that is to eigenstates which have one 
site occupied by two atoms and an extra empty site. Accord- 
ingly, these peaks are centered around the value hut = U. The 
reduced intensity of the 1-ph peaks with respect to the low- 
est energy peaks is expected, since in first-order perturbation 
theory the coupling of the ground state to 1-ph is proportional 
to 1/7. In fact as analytically calculated in Ref. |16||, for a 
unit filled system the height of the particle hole peaks is pro- 
portional to 64/7^ which is of the order observed in the plot. 
Panel (c) shows that the position and the height of the low- 
est energy peaks are reasonably well reproduced by the EF 
model, which here corresponds to standard fermionization, 
since n = 1. This is remarkable, since 7 is only 10 in Fig. 1(a), 
while it is considered infinite in the model. On the other hand, 
the high energy peaks are missing, since 1-ph are not present 
in the model llll . 

The qualitative analysis above is easily adapted to the case 
7\/ = 6, A^ = 7 of panels (b) and (d). Here, n = 2 and the 
EF model predicts that the low energy spectrum is the one of 
a single particle with hopping energy 2 J on top of a plateau 
of 6 atoms, which are frozen in a unit-filled Mott insulator. 
Thus, S{q,u!) in the EF model reduces to computing Eq. (|2ji 
for a non-interacting fermion with hopping energy 2 J. Since 
the single particle band width is now 8 J, we expect a large 
low frequency response in the energy range hui/ J < 8. This 
is actually shown in Figs, ^^b) and (d), which confirms the 
validity of the EF model even for such a small value of 7. In 
fact, the latter is just barely larger than 27c = 8. In addition, 
the exact results of Fig.fflb) show excitations around hui k, U . 
Not surprisingly, the latter are due to coupling to eigenstates 
with two sites occupied by two atoms, and to eigenstates with 
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FIG. 1: The dynamic structure factor S{q,Lu) vs u, with qd = 7r/3 
(continuous line), 27r/3 (dashed line) and tt (dotted line). Here, M = 
6 in all panels, while N — 5 (panels (a) and (c)) and N — 7 (panels 
(b) and (d)). Panels (a) and (b) are results of exact diagonalizations 
with 7 = U/J — 10, while panels (c) and (d) are the EF model 
results. For all plots, Q/J = 0. 
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a site occupied by three atoms. These couplings are again 
proportional to I/7, and thus the height of these high energy 
peaks is highly suppressed. We have also found numerically 
excitations at an energy of order 2U . Because the couplings of 
these eigenstates to the ground state are of order 1/7^, these 
peaks are very small, and are not shown in the graphs. 

Finally, we checked numerically that in the case M = 
N = 6, where the ground state is a Mott insulator with one 
atom per site, the lowest energy peaks are centered around U. 
The spectrum is thus gapped, in this case, as expected. The 
peaks' height is of the same order of magnitude of the peaks 
of Figs. 1(a) and (b) centered around hui JJ. A comprehen- 
sive discussion of the excitations in the MI phase can be found 
m Ref. d. In the following we show how the above discus- 
sion is useful in the description of the spectrum of excitations 
when the quadratic trap is present. 



V. INHOMOGENEOUS LATTICE (Q. > 0) 

When the quadratic trap is present and 7 > 7c, the system 
can show the coexistence of superfluid and insulating phases. 
Figure 12] displays the local density profiles of systems with 
such coexisting SF and MI phases, with onsite peak density 
one, (a), and larger than one, (b). In the plots, integer and non- 
integer local densities signal MI and SF phases, respectively. 
Here, AI = 100 and O = 0.008 J. In Fig.gja) the value of 7 is 
14, while it is 7 = 8 in Fig.|2jb). These values and the fillings 
we use are comparable to those used in current experiments. 
The continuous (red) and dashed (black) lines are the quan- 
tum Monte Carlo data and the predictions of the EF model, 
respectively. The two curves are almost indistinguishable on 
the scale of the graph, which is remarkable at least for the re- 
sults of panel (b), where the EF model at the center of the trap 
is barely applicable. For the case of Fig.|2ja) where the peak 
on-site density is one, the EF model corresponds to standard 
fermionization. On the other hand, for the case of Fig.lJJb), 
the EF model amounts to decomposing the density profile into 
two vertically stacked horizontal layers of atoms 1I14I1 (with 
A'^i = 61 and A^2 = 3 atoms in the lower and upper layers, 
respectively). In Fig.|2jb) atoms in the lowest layer have been 
shaded, in order to visualize the two layers of the EF model. 
In each layer there is at most one atom per site, and thus the 
EF idea is to apply standard fermionization techniques to each 
layer separately. Then, it is possible to calculate many body 
observables (e.g., the density profile) for the two layers sepa- 
rately, and add the results to obtain the total value. Details of 
the computation for the density profile are given in Ref. fl4ll . 
There, it is shown that EF provides accurate results whenever 
the superfluid regions at the trap center and at the edges are 
spatially well separated by sufficiently large MI regions. In 
fact, in this case it is possible to unanbiguously attribute a cer- 
tain number of atoms to each layer of the EF model, so that 
layers are well defined. This is the case of Fig.|2jb), where an 
insulating region of about 10 lattice sites separates the super- 
fluids at the trap center and at the edges. In the following we 
compute S{q,u!) for the two layers separately and then add 
the obtained results to calculate the complete S{q, uj). 
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FIG. 2: (Color online) Local density profile, rij, as a function of the 
lattice site j for the confined system with M — 100, and Q,/ J = 
0.008. Panel (a): iV = 50 and 7 = 14. Panel (b): iV = 64 and 7 = 
8. The continuous (red) and dashed (black) lines are the quantum 
Monte Carlo and EF model results, respectively. The shaded area in 
panel (b) corresponds to atoms in the lowest layer of the EF model. 
These atoms are largely frozen in a unit-filled Motts state, analogous 
to the case of panel (a). 
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FIG. 3: (Color online) ^(g, ui) vs ui for the system in Fig.|3a) with 
iV = 50, 7 = 14, and for qd = 0.27r(a), 0.57r(b), and 0.97r(c). 
The countinous(red) and dashed(black) lines are the quantum Monte 
Carlo and EF model results, respectively. The large excitation peaks 
for fttt;/ J < 4 are excitations of the superfluid component, while the 
broad peak at huj k, U = 14 J is due to particle-hole excitations of 
the combined SF and MI, analogous to Fig. 1 . 



Figure |3] shows S{q,Lo) as a function of oj for some val- 
ues of q for the system of Fig. 2(a), where there is at most one 
atom per site. In particular, Figs.|3ja), (b) and (c) correspond 
to qd = 0.27r, O.Stt and O.Ott, respectively. In the calculations 
we have normalized the height of the largest peak to 1 . Anal- 
ogous to the case of Fig. Q^a), the numerical solution shows 
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large peaks in the energy range hio < 4 J. The EF solution, 
which here reduces to standard fermionization, captures the 
location and width of these peaks rather well. Since fermion- 
ization accounts only for couplings to states with at most one 
atom per site, we conclude that the response for hw/ J < A 
is solely due to the small number of atoms in sites \j \ > 20 
which are in the superfluid phase (analogous to the discussion 
of Fig. ^a)). These are thus "surface" excitations. In addition, 
for q large enough {qd ~ tt) the numerical solution shows the 
presence of peaks around U, which are not captured by the 
EF model. These are due to couplings to 1-ph excitations, 
in analogy to the discussion above. Therefore, the excitation 
spectrum exhibits two branches, a lower energy one corre- 
sponding to excitations of the surface superfluid layer and a 
higher one corresponding to the Mott region. 

The differences between the numerical and EF results can 
be attributed both to errors of order ^ 1/7^ in the EF predic- 
tions and to the fact that the numerical calculation of S{q, uj) 
is difficult for 7 1. In this respect, we notice that the lowest 
energy excitation in the EF model is always at huj ^ UN ~ 
0.4 J, which is the energy cost for an atom at position N/2 
(farthest outlying occupied site), to tunnel to the next unoccu- 
pied site. The fact that the numerical results do not grab this 
excitation for qd ^ tt (Fig.|3lc)) suggests a partial failure of 
the numerics. 

Similarly, Fig.l^shows S{q, cu) (normalized to 1) as a func- 
tion of u! at different values of q for the case of Fig. |2jb), 
where there is a superfluid phase at the trap center with more 
than one atom per site. The continuous (red) line is the Monte 
Carlo result, while the dashed (black) line is the EF model 
result. As mentioned above, the EF model assumes that the 
system of Fig.lJJb) is conveniently decomposed into two dis- 
joint, vertically stacked gases, each one with at most one atom 
per site. Most atoms in the lowest layer are frozen in a Mott 
insulator state, in analogy to the case of FiglSja), while atoms 
in the upper layer are delocalized in a region approximately 
comprising the sites — 10 < j < 10. Atoms in the lower and 
upper layers have hopping energies J and 2 J, respectively, 
and are treated as independent. The results for S{q, ui) for the 
lower and upper layers are indicated by dotted lines (green 
and blue lines respectively). 

Figure|4]shows that the numerical excitation spectrum (con- 
tinuous line) at large enough momenta has two branches of 
excitations, which saturate at « 4 J and 8 J, respectively. 
The existence of these two branches of excitation has been 
previously reported by one of us in Ref. (rf\, where the lower 
and higher frequency peaks were attributed to excitations of 
the SF and MI phases respectively, as in the case of Fig. |3 
However, while this interpretation is correct for Fig. |3] it is 
only partially so in this case. The good agreement between 
the numerical results and the results of the EF model (dashed 
line) with respect to the position and width of the large fea- 
ture around huj « 4 J, and the position of the peak around 
huj « 8 J for large q, strongly suggests that this double peak 
structure is mainly due to excitations of the superfluid com- 
ponents of the two layers of the EF model. In particular, the 
lower frequency excitations {hui < 4 J) are due to superfluid 
atoms in the lowest layer of the EF model (sites where rij < 1, 
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FIG. 4: (Color online) S{q, lo) vs uj for the system in Fig.|2fb) with 
TV = 64, 7 = 8, and for qd = 0.27r(a), 0.47r(b), 0.77r(c), and 
0.97r(d). The countinous (red) and dashed (black) lines are the quan- 
tum Monte Carlo and EF model results, respectively. The dotted 
lines are results for the lower (green) and upper (blue) layers of the 
EF model. The dashed line is thus obtained by adding these curves. 
The excitation peak saturating at ficj/ J ~ 4 is due to excitations in 
the SF with less than one atom per site, (uj < 1, with 25 < \j\ < 35 
in Fig.|2lb)). The excitation peak saturating ai huj/J ~ 8 is due to 
excitations in the SF with more than one atom per site, (nj > 1, with 
< b l < 10 in Fig.|2lb))). Because 7 — 8, particle-hole excitations 
around ftt^j/J ~ 8 increase the width of this SF peak in the Monte 
Carlo results. 
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in the shaded region of Fig.|2ji, while the higher frequency ex- 
citations are mostly due to excitations of atoms in the upper 
layer, which are all in the SF phase. As mentioned above, the 
saturation thresholds 4J and 8 J correspond to the values of 
the effective single particle band-widths of atoms with hop- 
ping energy J and 2 J, respectively. 

Since here 7 = 8, 1-ph excitations certainly do contribute to 
the system's response around hcu/J^S. However, since the 
coupling to 1-ph excitations in the linear regime is rather sup- 
pressed, as can be seen in Fig.|3(it is proportional to 1 /j), we 
argue that the intensity of the response peak centered around 
hoj sa 8 J is dominated by the excitations of SF atoms in the 
second layer rather than by 1-ph excitations. Unfortunately, 
a detailed discussion of Fig.l^is complicated by the fact that 
the discrepancies between the model and Monte-Carlo results 
are most likely of physical as well as computational nature. 
As said above, the physical origin of the discrepancies resides 
in the EF model neglecting 1-ph excitations. This contributes 
to the underestimation of the height and width of the peak 
around ftw « 8 J in the EF results of Fig.^c) (while the peak 
around huj w 4 J is well reproduced). On the other hand, be- 
cause the relative height of the peaks in the numerical solution 
varies greatly from qd ~ O.Ttt (panel (c)) to qd = O.Ott (panel 
(d)) and oscillates in between (not shown), we are also led to 
suppose the existence of some instability in the numerical re- 
sults. This prevents us from attributing unambiguously all the 
discrepancies between the numerical and model results to 1- 
ph excitations, and thus getting a quantitative estimate of 1-ph 
contributions to the excitation spectrum. 

With all the provisos above, we can safely state that an 
important effect of coupling to 1 -ph excitations is to increase 
the width of the peak around fio; w 8 J with respect to the EF 
results. That is, by increasing the ratio 7, we would expect to 
see a continuous decrease of the width of the peak athoj/J = 
8, and the growth of a distinct peak at the chosen 7-value. 
This new distinct peak would be due to 1-ph excitations of 
the atoms in the two layers, analogous to the discussion of 
Fig.Ea-b). 

In other words, the presence of a multiply peaked structure 
in the low frequency response to the Bragg perturbation for 
7 ^ 71 J is clear evidence of the existence of a shell struc- 
ture in the many body density profile (Fig. |3- In fact, for 
qd ~ 11 the peak corresponding to the superfluid component 
of a generic layer n saturates at fio; w 4n J, and thus the pres- 
ence of low-frequency peaks separated by 4 J is an indica- 
tion of the cloud's shell structure. We thus propose to use 
Bragg spectroscopy to characterize the various phases of the 
trapped strongly correlated bosons, by detecting this multi- 
ple peak structure. For example, the system's response could 
be measured in successive experiments while keeping fixed 
the trapping potentials (and thus the 7 and Vl/ J ratios), and 



varying the total number of atoms N . Then, for « tt and 
a single peak at « 4 J would be observed 
{ ^/4U/n is a threshold value for having only one layer in 
the EF model, fl^ V By increasing N to values larger than 
^/AUjTl a peak at Sw/ J ?» 8 would be suddenly observed in 
the system's response, signaling the formation of a SF with 
two atoms per site at the trap center By further increasing 
N, other peaks would occur at multiples of 4 J, signaling the 
formation of superfluid regions with more than two atoms per 
site. It should be noted that typically the lifetime of these 
states would be severely limited by three-body recombination, 
and the observation of these layers may become complicated. 

Finally, we notice that Bragg spectroscopy in the linear 
regime may also be used to detect the atoms' shell structure 
in dimensions D larger than one (where fermionization tech- 
niques are not applicable). In fact, in the appropriate param- 
eter regimes, the system retains the structure of an array of 
vertically stacked horizontal layers. Thus, the response to 
the Bragg perturbation of atoms in the superfluid regions of 
the different layers should provide for a signature of the shell 
structure. Since the kinetic energy is a factor of D larger, 
the various branches of excitations at low frequency should 
be more separated for D > 1 than for D = 1. This may fa- 
cilitate the experimental resolution of the excitation peaks in 
dimensions higher than one. 



VI. SUMMARY 

We have presented a study of the dynamical structure 
factor for strongly interacting bosons confined in homo- 
geneous lattices and in presence of a quadratic confining 
potential. By comparing the results of quantum Monte Carlo 
simulations and exact diagonalizations with the predictions 
of an analytical EF model, we characterized the excitation 
spectrum of strongly correlated bosons at arbitrary densities. 
In particular, we were able to identify low frequency peaks 
in the system's response to the Bragg perturbation which are 
clear cut evidence of the formation of a shell structure in the 
density profile. This shell structure is made of alternating 
SF and MI phases, which we propose to study using Bragg 
spectroscopy in the linear response regime. 
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